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In order to improve aerodynamic performance and efficiency, monitoring the geometry and position informa- 
tion of the internal structure with the aero-engine under various operational statuses is an important task. In this 
paper, a novel design of industrial computed tomography (ICT) system with a linac as X-ray source is proposed 
to complete the task. The major advantage of the proposed system is that it can provide visualized images of 
internal structures of the running aero-engine without physical disturbance, which makes it possible to extract 
the accurate geometry information. The main idea behind this design is to measure the projection data from 
various views for reconstructions making use of the rotations of the aero-engine blades, instead of the mechan- 
ical rotations in the conventional ICT system. However, due to high speed rotation of the aero-engine blades, 
the system faces more challenges than conventional ICT systems both in data acquisitions and reconstruction 
algorithms. The challenges and corresponding solutions are presented in this paper. In conclusion, the proposed 
ICT system provides a powerful tool for monitoring the running aero-engines. 
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I. INTRODUCTION 


In aero-engine structure design and fabrication, running 
aero-engine test is an important experiment to improve the 
efficiency, aerodynamic performance and aircraft safety. For 
example, tip clearance has been a technical challenge in 
the aero-engine development. A 1% increase in tip clear- 
ance of the compressor casing will decrease the compressor 
efficiency by 1-3% [1], while probability of friction or col- 
lision between blade and casing increases with a too smal- 
l tip clearance, causing eventually serious accidents. Also, 
the blade dimensions affect greatly the engine aerodynamic 
performance. So, it is of great significance to monitor in- 
ternal geometry information of an aero-engine under various 
operation conditions. However, it is a harsh task to measure 
the inner structure under high temperature, high pressure and 
high rotating speed (10 000 rpm or higher) with conventional 
methods. For example, measurement accuracy of optical sen- 
sors can be reduced by oil dirt contamination [2], and eddy 
current measurement is restricted by the high temperature of 
aero-engines [3]. Also, it is not convenient or permissible to 
set up detectors or sensors in the aero-engines because they 
will disturb the airflow distributions. 

Currently, changes in dimensions of the blades under var- 
ious operating conditions are checked by indirect measure- 
ments, which cannot give a visualized image; or computer 
simulation, which cannot accurately simulate the complicat- 
ed real working environment. X-ray computed tomography 
provides accurate images without physical disturbance to the 
object being inspected, and is suitable to monitor the running 
aero-engine blades [4]. Another inherent advantage of X-ray 
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process tomography is hard-field character of X-ray, which 
makes it possible to achieve clear images without blurs or 
artifacts. Hampel et al. developed a gamma-ray tomogra- 
phy system to measure fluid distributions of a rapid rotat- 
ing fluid coupling [5]. They designed an angle synchronized 
trigger device to sort the projection data acquired from their 
system [6]. They succeed in revealing fluid structure of the 
coupling under different rotating speeds. However, an aero- 
engine is much more massive and complicated than that fluid 
coupling, and the numerous pipelines in the stationary part- 
s shall introduce significant contaminations to projections of 
the blades in high speed rotation. Also, determinations of the 
projection angles and suppressions of the motion blurs are 
challenges to the X-ray tomography because the blades rotate 
much faster than the fluid coupling. 

In this paper, an innovative design of industrial computed 
tomography (ICT) system is proposed for monitoring changes 
in internal structure of a running aero-engine. It shall make 
it possible to obtain accurate images of the internal structures 
under different operational statuses of the aero-engine. It has 
the potential to derive tomographic information about the w- 
hole engine running test. In Section II, details about the sys- 
tem design are presented. Three major challenges to the ICT 
system are analyzed and feasible solutions are proposed in 
Section III. In Section IV, we demonstrate feasibility of the 
ICT system and effectiveness of the solutions through com- 
puter simulations and experiments. Finally, a conclusion is 
given for the proposed ICT system and the major challenges. 


ll. METHODS 
A. System design 


Different from a static scanning object in a conventional 
ICT system, the aero-engine blades are rotating rapidly in the 
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process of measurements. An ICT system is developed to ac- 
quire images of rotating blades under steady operating state 
of an aero-engine. The main idea is to measure the projec- 
tion data by using the blade rotation of an aero-engine with 
stationary X-ray source and detector array [7, 8], rather than 
mechanical rotations of a conventional ICT system. As shown 
in Fig. 1, a 4MeV electron linac X-ray source produces a se- 
quence of X-ray beams pulsed at ~ 4 us width and the blades 
rotate at ~ 10000rpm. On the other side, the detector array 
receives the X-rays attenuated by the aero-engine. The data 
are transmitted to the workstation for image reconstruction. 
A sensor is used to monitor fluctuation of the blade rotation 
and the signals are used for accelerator control. 
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Fig. 1. (Color online) Schematic diagram of the proposed ICT sys- 
tem. 


In conventional ICT systems, the projections from differ- 
ent views are acquired in step-by-step mechanical turntable 
rotates of the CT system. The linac is triggered according to 
signal of angle encoders. However, in this design, it is the 
object that rotates. Two trigger modes (Fig. 2), single-trigger 
and multi-trigger modes, are proposed. In the single-trigger 
mode, the linac generates a sequence of X-ray pulses with an 
interval of to in several rotation cycles. The signals from the 
sensor are only used for estimating the rotation speed, but this 
is not necessary. This mode is more useful when the sensor 
signal is inexact or it is not permissible to set up such a sen- 
sor around the blades, because the trigger signal can be also 
sent manually. However, accuracy of the projection views is 
affected by the speed fluctuation in this mode. In the multi- 
trigger mode, the sensor sends a pulse signal to the system 
controller whenever a blade passes by. The accelerator is syn- 
chronized to the blade rotations. Particularly, the controller 
make the linac triggered with a proper delay of t,, and a peri- 
od of Tperioa according to the rotation speed and ideal projec- 
tion angles. A programmable logical controller (PLC) is used 
to achieve this function. In this mode, uniform samplings are 
available even though the rotation speed is not constant, be- 
cause the controller adjusts the delay time t„ according to the 
speed fluctuation. 


B. Reconstruction algorithm 


The dynamic measurement procedure and the rotating 
blades make the recovery of the internal structure more com- 
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Fig. 2. (Color online) The sequence diagram of the two trigger mod- 
e: single-trigger and multi-trigger mode (M = 4 pulses in each 
cycle). 


plex and difficult. There are three major obstacles to over- 
come. First, it is crucial to sort the projections with corre- 
sponding projection views and ensure quality of the recon- 
structed image. Next, the “blocked projections” must be cor- 
rected. As all projections are acquired by using the rotat- 
ing blades, the stationary casing and pipelines around the 
blades manifest as block contaminations in the projection da- 
ta, which brings significant artifacts in the final results. Fi- 
nally, there are the motion artifacts. Due to the high rotation 
speed, the blades rotate a certain angle in an X-ray exposure, 
hence the blurred edges. In the following sections, the three 
challenges will be analyzed. 


I. Determination of projection views 


In conventional ICT systems, the exact projection views 
can be achieved from the angle encoder set up on the mechan- 
ical turntable. However, in this design, the angle sampling in- 
terval is determined by the product of rotation speed and the 
time interval between two adjacent X-ray pulses. Determina- 
tion of the exact projection angles is the key to achieving high 
quality images. In the single-trigger mode, the linac produces 
pulsed X-ray beams in a given frequency. The sampling angle 
can be estimated according to the accelerator frequency and 
blade rotation speed if the latter is constant. However, due 
to fluctuation of the rotation speed, the projection angles are 
often slightly inaccurate and the errors accumulated in mea- 
surements of several rotation cycles would eventually lead to 
obvious stripe artifacts in the reconstructed images. In order 
to reduce the accumulated errors, a posterior correction-based 
method with statistical properties of the projections is pro- 
posed to solve the problem. The projections measured from 
the same view or two views quite close to each other should 
have a very high similarity. The similarity can be measured 
by cross-correlation coefficient or by standard deviation. For 
example, in the single-trigger mode, after M rotating cycles, 
the N" X-ray pulse is produced when the blades rotate to the 
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same position where the accelerator generates the 1%* X-ray 
pulse. Then the projection must be quite similar to the 1“ 
projection and the angle sampling interval 69 can be estimat- 
ed by Eq. (1) and the n" projection views can be calculated 


by Eq. (2). 
bo = 360M /N, (1) 


On = (n — 1)9o. (2) 


From Eq. (2), the error of projection views increases with 
n. In order to suppress the accumulated errors, adjacent pro- 
jections are grouped and the projections are divided into sev- 
eral groups. Then images are reconstructed with these projec- 
tions, respectively. While limited amount of projections sup- 
presses the accumulating errors, insufficient projection data 
may affect the image qualities. For taking full use of all the 
projections, registrations between these pre-reconstructed im- 
ages are done to acquire the relations of projection angles 
(offset angles) among different projection groups. Finally, all 
the projection data are reorganized for reconstructions while 
the accumulating errors are greatly limited. The complete 
work flow is presented in Fig. 3. 
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Fig. 3. (Color online) Work flow of the posterior correction based 
method. The scale bar is set to [0.0 1.0]. 


With a sensor to detect the blades, the multi-trigger mode 
can be adopted, and uniform angle samplings can be achieved 
by the controller, in which a PLC adjusts the delay time and 
pulse frequency adaptively to match the changing rotation 
speed. Projection views incorporating information from the 
sensor are more accurate. If the speed change is too big, the 
projections obtained far from the ideal positions can be re- 
moved or re-measured. 


2. Correction of blocked projections 


The block contaminations in the projections data are 
caused by the stationary casing and complicated pipelines 
around the blades. It is noted, however, that both of the 


Nucl. Sci. Tech. 25, 060202 (2014) 


(d) 


(a) Reconstruct the image using the contaminated projection; 


(b) An interaction-based registration between the prior image 
and the reconstructed image; 


(c) Estimate the block components and correct the projection 
by Eqs.(4) and (5), respectively; 


(d) Reconstruct the image of objects under operating test with 
the reconstructed data. 


Fig. 4. (Color online) Work flow of the interaction based method 
with a prior image. The scale bar is set to [0.0 0.04]. 


source-detector and block parts keep static in the whole 
scanning procedure, which means that the block components 
to a specific detector at different views are constant. The ex- 
pression of discrete imaging model at i projection view to 
the j detector can be given as 


Ay jf +b; = ij, (3) 


where H is the projection operator, f is the image datum, b 
is the block term, and p is the projection datum. Details of 
the reconstructed image using the contaminated projections 
are seriously corrupted by these artifacts, making it hard to 
determine the blade edges. 

An interaction-based method with a prior image is devel- 
oped to address the problem. Contributions of the block com- 
ponents from the acquired projection data will be removed, or 
they will be estimated. In this paper, a prior image is incorpo- 
rated to achieve the goal. Before the dynamic measurements, 
a Static aero-engine is scanned by a conventional ICT and an 
image without blocks is reconstructed. This can be used as 
the prior image fprior in the procedure. The prior image fprior 
is used to simulate the contribution of H; ; f for estimating 
the block components though there are some differences in 
details. It should be noted that the registration between the 
reconstructed image and the prior image is essential before 
the estimation, because they are generally not at the same lo- 
cation. Some serious artifacts besides the deformations can 
appear in the reconstructed image, making it difficult to apply 
auto match technologies, such as the Scale-Invariant Feature 
Transform (SIFT) method [9], for the registration. Hence, 
an interaction is introduced into the registration procedure. 
Then, a forward-projection of the matched prior image fprior 
with the practical scanning configuration parameters is per- 
formed. The block components of gh detector, bj, can be es- 
timated by 


= H; j fprior- (4) 
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By Eq. (5), the b; are used to obtain corrected projection data, 
which are used to reconstruct the final image. 


Diz = Pij — bj. (5) 


The summarized work flow is shown in Fig. 4. 


3. Suppression of exposure blurs 


The exposure blurs in the final reconstructed images are 
due to the speed high rotation, by which the blades rotate a 
certain angle during an X-ray exposure. This gives rise to 
significant degradations of the image spatial resolution and 
great limitations of the ICT in high-precision measurements. 
Therefore, the blade rotations during the X-ray exposure must 
be considered in the reconstruction method. An iterative re- 
construction method is proposed based on the discretization 
of the X-ray exposure procedure. The main idea behind the 
method is to discrete the dynamic X-ray exposure to several 
virtual conventional exposure child procedures. Specifically, 
this physical procedure is divided into M virtual convention- 
al static exposure processes. In different sub-processes, the 
blades rotate to different positions as 


1 M M 


= MI;/Io, 


where the subscripts j and i indicate the j} X-ray path and 
the number of sub-process of virtual conventional static ex- 
posure, respectively. In fact, H; j is the system matrix of ge 
detectors at i® sub-process. 

Equation (6) is a nonlinear problem raising great chal- 
lenges to the conventional reconstruction methods, where the 
scanned object is assumed static. Approximately, however, it 
can be converted to a linear problem using the Taylor expan- 
sion with a good initial guess fo. 


M 
H; = > exp (Huy fo) Hil, 
i=l (7) 


MI; 


M 
Pj = X exp (—H;; fo) = 
i=1 

We note that the difference image df instead of f is calcu- 
lated in each iteration. Therefore, the initial guess fp should 
be updated after each iteration by fo = fo + df. 

According to our research, the initial guesses can be from 
either a prior image or the blurred image. Then, it can be 
solved by widely used algebraic reconstruction methods on- 
ly with some simple modifications on the system matrix and 
projection data as shown in Eq. (7). 

In summary, Fig. 5 shows the total flowchart for a com- 
prehensive demonstration of the system challenges and the 
corresponding solutions. 
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Fig. 5. The total flowchart of the system challenges and correspond- 
ing solutions. 


Ill. SIMULATIONS AND EXPERIMENTS 
A. Numerical simulations 


In computer simulation, feasibility of the single-trigger 
mode of the proposed system was validated first. As shown 
in Fig. 6(a), a simple phantom is designed to simulate the 
aero-engine blades. The angle-intervals of samplings were 
estimated based on similarities of projections as mentioned 
above, but errors of the estimated angle-intervals were in- 
evitable. The errors accumulated in projection angles be- 
cause projection angles were calculated by accumulating the 
angle-intervals. As shown in Fig. 6(b), quality of the image 
reconstructed with rough projection angle positions degrad- 
ed, which greatly influenced the following quantitative anal- 
ysis. On the contrary, visible improvements in image quality 
were achieved using the proposed posterior correction based 
method. The relative error reduced to 48% of the convention- 
al one because the maximum estimated errors of projection 
angles decreased to 1.61° from 4.3°. Hence, the posterior 
correction of projection angles is essential to guarantee the 
image quality in single-trigger mode when the proposed sys- 
tem is applied to monitor the running aero-engine blades. 


Fig. 6. Simulation results with simple phantom. (a) Original im- 
age; (b) image reconstructed with rough angle positions; (c) image 
reconstructed by posterior correction based method. The scale bar is 
set to [0.0 1.0]. 
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Next, an air blower was used as a phantom to demonstrate 
the effectiveness of the method for suppressing the exposure 
blurs. However, air blowers rotate at about 2900 rpm, rather 
than 10 000 rpm of aero-engine blades. In order to simulate 
projections of real rotating aero-engine blades, the successive 
projections were accumulated to simulate significant rotating 
movements during the X-ray exposures. The results of con- 
ventional method and ICT are presented in Fig. 7. 


Fig. 7. Simulation with an air blower as the phantom. (a) Reference 
image; (b) image reconstructed by conventional method; (c) and (d) 
images reconstructed by our method with 8 and 16 virtual exposures, 
respectively. The scale bar is set to [0, 0.04]. 


In Fig. 7, compared with the reference image, significan- 
t exposure blurs and the decrease of contrast-to-noise ratio 
(CNR) can be observed in Fig. 7(b) using conventional re- 
construction method due to the rotations during the X-ray 
exposure, while the images reconstructed by our method are 
of sharp edges (Figs. 7(c) and 7(d)). For quantitative anal- 
ysis, the CNR at blade edges were calculated. The CNRs 
of Figs. 7(b), 7(c) and 7(d) were 50, 70 and 80, respective- 
ly. So, the more virtual child exposures used to simulate the 
actual exposure, the better result will be achieved, with heav- 
ier computations, though. In practice, the number of virtual 
exposures is often determined by the tradeoff between image 
quality and computations. The development of corresponding 
accelerated method using GPU will be part of future work be- 
cause of its high parallelism. The high quality image makes 
it possible to detect minor changes in internal structure of the 
running aero-engine. It also demonstrates the effectiveness of 
our method in suppressing exposure blurs. 


B. Experiments 


In order to demonstrate feasibility of the proposed ICT sys- 
tem and effectiveness of the methods, a 2 900 rpm air blower 
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Fig. 8. (Color online) The air blower used in the experiment simu- 
lating the aero-engine. 


was scanned. Seated on the casing, a laser sensor generat- 
ed a pulse when the blade tip passed by. In order to simu- 
late pipelines of an aero-engine, several steel columns were 
placed randomly around the air blower (Fig. 8). In this ex- 
periment, multi-trigger mode was adopted. The maximum 
frequency of the 4 MeV linac is 250 Hz, so 4 samplings were 
acquired in each cycle. The total 1 024 measurements were 
completed in 256 cycles. In other words, the whole scanning 
time is in less than 6s, which is significantly faster than con- 
ventional CT measurements. It is possible to achieve several 
tomographic images in the test. 

Firstly, a conventional ICT scan with the air blower in 
static was performed and internal structure of the air blow- 
er under idle state was reconstructed (Fig. 9(c)). Next, dy- 
namic measurements were carried out with the air blower 
rotating at about 2 900 rpm. The internal structure was re- 
constructed with the reorganized projection data (Fig. 9(b)). 
However, some significant ring artifacts caused by the steel 
columns corrupted the image quality (Fig. 9(d)). Then, the 
projection data were corrected by using the interaction-based 
method with the static image as the prior image. It can be 
seen from Figs. 9(e) and 9(f) that the reconstructed image of 
blade shapes differs significantly from the prior image. So, a 
registration was carried out according to the major structures, 
such as the air blower casing. The final image (Fig. 9(g)) re- 


060202-5 


CHANG Ming, XIAO Yong-Shun and CHEN Zhi-Qiang 


Nucl. Sci. Tech. 25, 060202 (2014) 


Fig. 9. (Color online) Experiment results with an air blower. (a) original sinogram; (b) reorganized sinogram; (c) reconstruction of conven- 
tional static ICT; (d) image reconstructed using the projections without correction; images ((c) in red channel and (d) in blue channel) before 
(e) and after (f) registration; (g) image reconstructed using the projections with correction. The scale bar is set to [0, 0.04]. 


constructed using the projections corrected by the interaction- 
based method was almost free of ring artifacts, and the main 
differences in blade shapes were well recovered. These im- 
provements in the image quality are meaningful for geomet- 
ric measurements. But there were still some local artifacts 
caused by the systematic errors, such as fluctuation of the sen- 
sor signal, mechanical vibration and some other engineering 
problems. The experiment design with a sensor of enhanced 
stability will be adjusted for further improvement of the im- 
age quality. 


IV. CONCLUSION 


In this paper, an innovative design of X-ray process to- 
mography system is presented to obtain internal structure im- 
ages of running aero-engine blades, which makes it possible 
to extract accurate geometry information of internal struc- 
tures of an aero-engine under different operational statuses. 
Compared with conventional measurements, it has the po- 
tential to achieve high-quality image under harsh measure- 


ments without physical disturbance. The image sequence 
about each steady state during the running aero-engine test 
can be achieved by the proposed system, which is meaning- 
ful to the design and control of the aero-engine, though only 
certain steady state is involved in this paper. Its advantages 
of non-contact and visual imaging make it suitable for moni- 
toring the internal structure changes of running aero-engines. 
Three major challenges caused by the dynamic measurements 
and rapid rotating movements are discussed. Some feasible 
solutions are raised in this paper and the effectiveness of these 
methods is demonstrated by simulating and experimental re- 
sults. These basic researches are meaningful to the following 
system implementations in the future. Although the prob- 
lem of measuring tip clearance was not solved completely, 
the proposed system with corresponding algorithm provides 
anew powerful tool based on CT technology for not only dy- 
namic measurements of the tip clearance but also inspections 
of the internal geometry and position changes. Furthermore 
our system design and solutions are not limited to monitor 
running aero-engine blades. It can be also further expanded 
to other X-ray process tomography applications in the fields 
of surveying internal structure changes of rapid moving ob- 
jects, such as gas turbines or steam turbines. 
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